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Two, non-interacting systems immersed in a common bath and evolving with a Marko- 
vian, completely positive dynamics can become initially entangled via a purely noisy 
mechanism. Remarkably, for certain, phenomenologically relevant environments, the 
quantum correlations can persist even in the asymptotic long-time regime. 

1. Introduction 

The simplest way to generate quantum correlations between two systems is via 
a suitable Hamiltonian coupling: optimization of entanglement production can be 
achieved by a careful choice of the form of the interaction term. 1-5 

When the two systems are immersed in an external bath, decoherence phenom- 
ena usually occur, counteracting entanglement generation. These effects are clearly 
a curse in quantum information and in fact various error correcting strategies have 
been devised in order to limit these environment induced mixing enhancing phe- 
nomena. 

However, an external environment can also provide a further, indirect cou- 
pling between the two systems and therefore an additional mechanism to correlate 
them. 6-9 That the external environment can indeed generate entanglement has 
been first established in exactly solvable models: 6 there, correlations between the 
two subsystems take place during a short time transient phase, where the reduced 
dynamics of the subsystems contains memory effects. 

Remarkably, a similar phenomenon of entanglement production may occur also 
in the Markovian regime, through a purely noisy mechanism. 10-12,28 ' 29 For the 
case of two, two-level systems immersed in a common bath this phenomenon can be 
established by looking at the eigenvalues of the partial transposed density matrix 
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that represents the two subsystem state. 13,14 A sufficient condition for initial, envi- 
ronment induced entanglement generation can then be obtained: it allows deriving 
a test on the entanglement power of the bath. 10 

Nevertheless, this test is unable to determine the fate of the initially created 
quantum correlations as time becomes large. In order to discuss asymptotic entan- 
glement, one has to analyze directly the structure of the dissipative, Markovian 
dynamics followed by the two systems, and determine its equilibrium states. 

In the following, we shall present such an investigation for a subsystem com- 
posed by a couple of identical two-level systems, evolving with a completely positive 
dynamical semigroup. Being interested in discussing the correlation power of the 
environment, we shall assume the two systems to be independent, without any mu- 
tual direct interaction^ As we shall see, there exists a class of environments, which 
is of relevance in phenomenological applications, that are able not only to initially 
generate entanglement: they can continue to enhance it even in the asymptotic long 
time regime. 



2. Master Equation 

We shall deal with a subsystem composed by two, identical, non-interacting two- 
level systems, immersed in a common heat bath. The derivation of a physically 
consistent master equation for the reduced density matrix p(t) = Tr^ptot^L °b~ 
taincd by tracing the total density matrix p to t (t) over the bath degrees of freedom, 
is notoriously tricky, requiring an a priori unambiguous separation between subsys- 
tem and environment. 15-18 ' 28 Generally speaking, this distinction can be achieved 
when the correlations in the environment decay much faster than the characteristic 
evolution time of the subsystem alone. Then, in the limit of weak couplings, the 
changes in the evolution of the subsystem occur on time scales that are very long, 
so large that the details of the internal environment dynamics result irrelevant. 

This situation is amenable to a precise mathematical treatment: 19 as a result, 
the two-system state p(t) evolves in time according to a quantum dynamical semi- 
group of completely positive maps, generated by a master equation in Kossakowski- 
Lindblad form: 19,20 ' 21 

^ = -i[n eS ,p(t)]+jr[p(t)}. (1) 

The unitary term depends on an effective Hamiltonian, containing both the initial 
system Hamiltonian and suitable Lamb contributions; in general, it can be decom- 
posed as: Heft — #eff + -^iff + Heft * 1 ■ The first two terms represent single system 
contributions, 

3 3 

ff^E^fVi®!), ffW= ^ ^(KSaO, (2) 



i=i 



a In the same way, we shall also neglect any Lamb shift Hamiltonian contribution that might be 
generated by the presence of the external bath. 



Asymptotic Entanglement of Two Independent Systems in a Common Bath 3 



where (<7j (g) 1), (1 <g> (Tj) are the basis operators pertaining to the two systems, 
respectively, with Uj, i = 1, 2, 3 the Pauli matrices; the third piece is a bath generated 
direct two-system coupling term, that can be expressed as: 

^ = XX^fo®^-). (3) 

The dissipative contribution £[/»(£)] can be cast in Kossakowski form, 20 

6 



Q,/3=l 



(4) 



using the hermitian, tracelcss, matrices .F Q that coincide with the first system basis 
operators (a a ® 1) for a = 1, 2, 3, while reproducing the second system basis opera- 
tors (1 ®<T a _3) for a = 4, 5, 6. The Kossakowski matrix C a p is a 6 x 6 matrix, which 
is non-negative, thus guaranteeing the complete positivity of the reduced dynamics. 
Using the splitting of the indices introduced above, it can be conveniently written 
as 

c-(iVi, OS 



& C , 

in terms of 3 x 3 matrices A = A\ C = and B. This decomposition carries 
a direct physical interpretation. Indeed, the pieces in (4) containing the diagonal 
contributions A and C correspond to noise terms that affect the first, respectively 
the second, system in absence of the other. On the contrary, the pieces depending on 
B encode environment generated dissipative couplings between the two, otherwise 
independent, systems. 

In order to obtain a more explicit expression for the dynamics, it is convenient 
to decompose the 4x4 density matrix p{t) along the Pauli matrices: 



Pit) = \ 



(6) 



where the coefficients poi (t) , p^ (t) , p^ (t) are all real. Substitution of this expansion 
in the master equation (1) allows deriving the corresponding evolution equations for 
the above components of p(t) . As mentioned in the Introduction, we are interested in 
studying possible entanglement production through the purely dissipative action of 
the environment; in the following, we shall therefore ignore the Hamiltonian pieces 
and concentrate on the study of the effects induced by the dissipative part C[ ■ ] in 
(4) b 



b In other terms, we limit our analysis to baths for which the induced two-system Hamiltonian 
coupling in (3) is vanishingly small or alternatively does not give rise to temperature dependent 
entanglement phenomena. This situation is rather common in phcnomcnological applications; for 
concrete examples, see Refs.[ll,12]. 
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Although the general case can be similarly treated, for sake of simplicity we 
shall here limit our considerations to baths for which the submatrices in (5) are 
all equal: A = B = C. This choice of the Kossakowski matrix, although special, is 
nevertheless of great phenomenological relevance since it is adopted in the analy- 
sis of the phenomenon of resonance fluorescence. 22 ' 23 In addition, it is precisely a 
dissipative term of this type that describes the interaction of two atoms with a set 
of weakly coupled external quantum fields, in the limit of a vanishing spatial atom 
separation. 11 ' 12 

In such a situation, the form of the dissipative contribution in (4) simplifies so 
that the evolution equation can be rewritten as 



dp(t) 
dt 



p{t) 



1 



Pit)} 



in terms of the following symmetrized two-system operators 
S 4 = a, 1 + 1 ®g % , i = 1,2,3 . 



(7) 



(8) 



One easily checks that these operators obey the same su(2) Lie algebra of the Pauli 
matrices; further, together with 



' a 3 + a 3 ® a i i h 3 — 1) 2 , 3 , 



(9) 



J 3 ' "3 

they form a closed algebra under matrix multiplication, whose explicit expression 
is collected in the Appendix. 

It is now convenient to decompose the hermitian matrix Aij into its real and 
imaginary parts, 

3 

Aij = A^ + i ^2 SijkB k , (10) 
fe=i 

with A^ real symmetric and Bi real. Inserting this in (7) and using the decom- 
position (6), a straightforward but lengthy calculation allows derive the following 
evolution equations for the components of p{t) : 

3 

-2Ap 0i (t) + 2Y^ \AikPok{t) - Pik(t)B k ^ +2(2 + r)B 4 , 



dt 

dpiojt) 
dt 

dt 



k=l 
3 

-2Ap i0 (t) + 2Y^ [A ik pko{t) -pki(t)B k ] + 2(2 + T )B t , 
fe=i 

3 



AA[ Pij (t) + pji(t)] +"2^2 [A ik pkj(t) + A jk p ik (t) 



k=i 



AA i3 T 



<3 <3 

+4 5Z [ A ikPjk(t) + Aj k p ki (t) +4 At- A u pik{t) 



k=l 



k,l=l 



+2[Bip j0 (t) + Bjp i(t)j +A\B iPoj {t) + BjPio(t)j 
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3 

-2^B k [p ak (t)+p k0 {t)] Sij . (11) 

In these formulas, the parameter A represents the trace of while the quantity 
r = J2i=i pu that of the submatrix pij. By taking the trace of both sides of the 
last equation above, one discovers that r is a constant of motion. Nevertheless, the 
value of r can not be chosen arbitrarily; the requirement of positivity of the initial 
density matrix p(0) readily implies: — 3 < t < 1. 

3. Environment Induced Entanglement Generation 

Using the explicit form (11) for the master equation derived in the previous Section, 
one can now investigate whether an external environment can actually entangle the 
two independent systems. Since we are dealing with a couple of two-level systems, 
this can be achieved with the help of the partial transposition criterion: 13,14 a state 
p{t) results entangled at time t if and only if the operation of partial transposition 
does not preserve its positivity. 

We shall first discuss the possibility of entanglement creation at the beginning of 
the evolution: if a bath is not able to initially entangle the two systems, it will hardly 
do so in the limit of large times. A simple strategy to ascertain entanglement creation 
is as follows: assume the initial state to be pure and separable, i.e. p(0) = \<p)(<p\ ® 
\ip)(ip\, and then find out whether in the neighborhood of t = the dynamics in 
(11) is able to make negative an initially zero eigenvalue of the partially transposed 
density matrix p(t) (note that p(0) = p(0) for the chosen initial state). This amounts 
to study the behavior of the time derivative d t p(0), that can be explicitly obtained 
by taking the partial transposition of both sides of (7) (or equivalently of the system 
of equations in (11)). 

In this way, one finds that the dynamics generated by the equations in (7) can 
indeed entangle the two subsystems, provided at least one of the coefficients Bi 
in the Kossakowski matrix Aij in (10) is nonvanishing (see Ref.[10] for further 
details). This is a generic property of the Markovian dynamics in (7): entanglement 
is generated as soon as t > 0. 

In order to study the fate of this initially produced correlations as time becomes 
large, one needs to analyze the ergodic properties of the semigroup evolution gen- 
erated by (7). On general grounds, one expects that the effects of decoherence and 
dissipation that counteract entanglement production be dominant at large times, 
so that no entanglement is left at the end. However, as we shall explicitly see, there 
are cases for which the environment induced entanglement creation never stops as 
time flows, allowing at the end the presence of entangled equilibrium states. 



c By analyzing the structure of a general evolution equation with the dissipative term as in (4), 
one can show that this is the case only when the condition A = B = C is satisfied. This result is 
also related to the existence of multiple equilibrium states; see below. 



6 F. Benatti and R. Floreanini 



The system of first order differential equations in (11) naturally splits into two 
independent sets, involving the symmetric, p/ 0i ) = p i + Pm, Puj) = Pij + Pji, 
and antisymmetric, pr oi i = poi — p^, puj\ = pij — Pji, variables. By examining 
the structure of the two sets of differential equations, one can conclude that the 
antisymmetric variables involve exponentially decaying factors, so that they vanish 
for large times. Then, recalling the definitions in (8) and (9), the study of the 
equilibrium states p of the evolution (7) can be limited to density matrices of the 
form: 



Pij Sij 

i—l i->3 = 1 



(12) 



with p^ = pji. 

These states obey the equilibrium condition d t p = 0, and therefore annihilates 
the r.h.s. of all equations in (11). By direct inspection, one finds that these con- 
ditions are invariant under linear orthogonal transformations that act on both the 
coefficients Bi 1 Aij and the components pi, pij of the density matrix p. Then, with- 
out loss of generality, for the purpose of identifying the asymptotic states, one can 
take the real part of Aij to be diagonal, i. e. = Aj Sij ; the general case can always 
be recovered at the end by undoing the orthogonal transformation that has brought 
A^ in diagonal form. 

Further, in order to simplify the exposition, we shall assume the vector of com- 
ponents Bi to be directed along the third axis, so that only the component B 3 = B 
will be nonvanishing. Then, positivity of Aij readily implies: Aj > 0, i — 1,2,3 and 
B 2 < AiA 2 . 

The approach to equilibrium of semigroups whose generator is of the generic 
Kossakowski-Lindblad form has been studied in general and some rigorous math- 
ematical results are available. 24 ' 16 We shall present such results by adapting them 
to the case of the evolution generated by the equation (7). 

First of all, one notice that in the case of a finite dimensional Hilbert space, there 
always exists at least one stationary state pa. d Let us now introduce the operators 
Vi = A 1 / 2 ( recau tlmt A is non- negative) , so that the r.h.s. of (7) can be 

rewritten in the so-called diagonal form: 



3 

(13) 



When the set M formed by all operators that commute with the linear span of 
{Vi, Vj, i — 1,2,3} contains only the identity, one can show that the stationary 
state po is unique, and of maximal rank. On the other hand, when there are several 



d This can be understood by recalling that in finite dimensions the ergodic average of the action 
of a completely positive one-parameter semigroup on any initial state always exists: the result is 
clearly a stationary state. 
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stationary states, they can be generated in a canonical way from a po with maximal 
rank using the elements of the set M. 

In the present case, Ai contains the operator S = Su, besides the identity; 

indeed, with the help of the algebraic relations collected in the Appendix, one imme- 
diately finds: [S, Sj] = 0. Out of these two elements of M, one can now construct 
two mutually orthogonal projection operators: 



4 



2 



Q = l-P 



(14) 



Then, one can show that any given initial state p(0) will be mapped by the evolution 
(7) into the following equilibrium state: 



P(0) 



PpoP 
Tr[Pp P] 



Tr[Pp(0)] 



^[Qpo Q] 



Tt[Q P (0)] 



(15) 



That this state is indeed stationary can be easily proven by recalling that P and Q 
commute with Ej, i = 1,2, 3, and thus with the Vi as well; therefore, C[p] = 0, for 
any p(0), as a consequence of £[po] = 0. 

The problem of finding all invariant states of the dynamics (7) is then reduced 
to that of identifying a stationary state po with all eigenvalues nonzero. Although 
in principle this amounts in solving a linear algebraic equation, in practice it can 
be rather difficult for general master equations of the form (1). Nevertheless, in the 
case at hand, the problem can be explicitly solved, yielding: 



Po 



1<S>1 + M£ 3 -N(S U -S22) +RS; 



:;:! 



with 



M 



N 



R = 



2B 



Ai + A2 



(Ai -\ 2 )B 2 



2(Ai + A 2 )(A 1 A 2 + A 1 A 3 + A 2 A 3 ) ' 

(Ai + A 2 + 4 A 3 ) B 2 
2(Ai + A 2 )(AiA 2 + A 1 A 3 + A 2 A 3 ) ' 



(16) 

(17) 
(18) 
(19) 



Note that the previously mentioned positivity conditions on the parameters Ai, 
i = 1,2,3 and B, put restrictions on the components of po given above; in particular, 
one has: < 2 R < 1, M 2 < 2R, M 2 +4N 2 < 1, and the upper limits can be reached 
only when Ai = A 2 . f 

Inserting this result in the expression (15) allows deriving the expression of 
the set of all equilibrium states of the dynamics (7); as expected, they take the 



e One easily checks that P is the projection operator on one of the maximally entangled Bell states. 
f These inequalities also guarantee that po as given in (16) is indeed a state, i.e. that all its 
eigenvalues are non-negative. 
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symmetric form of (12), with the nonvanishing components given by: 

h = ^M, (20) 

(1 + 2N)t + 2(3N-R) 

911 = 2(3T^ ' (21) 

_ (l-2N)r-2(3N + R) 

922 2(3 + 2R) ' (22) 

933 ~ 2(3 + 2R) ■ (23) 

These stationary density matrices depend on the initial condition p(0) only through 
the value of the parameter r, that as already mentioned is a constant of motion for 
the dynamics in (7). 

Now that we have completely classified the stationary states, one can study their 
properties, in particular with respect to quantum correlations. It turns out that p 
is in general entangled. 

To explicitly show this, one can as before act with the operation of partial 
transposition on p to see whether negative eigenvalues are present. Alternatively, one 
can resort to one of the available entanglement measures and concurrence appears 
here to be the more appropriate: its value C[p] ranges from zero, for separable states, 
to one, for fully entangled states. 25-27 In the case of the state p above, one finds 



C[p] 

where 



JJ 



[2(3 + 2R) 



(1 - 2R) 2 + 4(2i? - M 2 ) "' " . (25) 



4R- 3A 
2 + A 



, (24) 



1/2 



The expression in (24) is indeed nonvanishing, provided we start with an initial 
state p(0) for which 

4-R-3A 

t < — — . (26 

2 + A y ' 

The concurrence depends linearly on the initial parameter r; it assumes its maxi- 
mum value C[p] — 1 when r = —3, as for the state P in (14), and reaches zero at 
t = (AR - 3A)/(2 + A) < 1. 

This result is remarkable, since it implies that the dynamics in (7) not only can 
initially generate entanglement: it can continue to enhance it even in the asymptotic 
long time regime. In other terms, prepare the two atoms in a separable state at t = 0; 
then, provided the condition (26) is satisfied, their long time equilibrium state will 
turn out to be entangled. 

Entanglement enhancement is not limited though to initially separable states: 
one can easily check that the phenomenon of entanglement production through a 
purely noise mechanism takes place also when the initial state p(0) already has a 
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non- vanishing concurrence. As an example, let us consider the following initial state, 
built out of the two projector operators introduced in (14): 

P(0) = |Q + (1 - s)P ; (27) 

it interpolates between the completely mixed (separable) state obtained for s = 3/4 
and the totally entangled state P. Provided s < 1/2, this state is entangled, with 
C[p(0)] = 1 — 2s. The difference in concurrence as this initial state evolves to its 
corresponding asymptotic one p turns out to be 



C[p]-C[p(0)} = 2s 



I- 2 



3 + 2R 



(28) 



which is indeed non vanishing. As a final remark, notice that this enhancement in 
concurrence vanishes as s approaches zero; in other terms, the maximally entangled 
state P can never be reached as an asymptotic state unless one already starts with 
it at t = 0: P results an isolated fixed point of the dynamics generated by (7). 

Appendix 

We collect here the algebraic relations obeyed by the nine hermitian, traceless ma- 
trices Sj = <7j (g> 1 + 1 ® &i and Sij = Gi <g> a j + Oj <S> a i , i,j — 1,2,3, introduced in (8) 
and (9). As mentioned in the text, altogether they form a closed algebra under ma- 
trix multiplication. In fact, using <Ji(jj = Sij + i^/s=i E ijk a k, a direct computation 
yields: 



Ej = 2 6jj 1 ® 1 + i e^k Sfc + >% , 

3 3 

Sij Sfe = 5 t k Xj + Sjk E, + i ^ Siki Sij + i ^ Eju Su , 

i=i i=i 

3 3 

Sfe Sij = 5ik Ej + Sjk Ej - i £ iki Si 3 - i ^ Ejki Su , 

i=i i=i 

3 

Sij Sfcj = 2((5j fe 5j/ + <5j; <g>l + i^2 ($ik£jlr + SjkSilr + SuSjkr + SjlSikr^J S r 

r=l 

-^2% 4; - ^ife (5j/ - Su Sjkj S + 2(s i3 S k i + S k i S i3 ^j 
—Sik Sji — Su Sjk — Sjk Su — Sji Sik , 

where S — ^2 r=i Su. From these relations, one immediately sees that the commutant 
of the linear span of the set {Ej : i, j = 1, 2, 3} contains two elements, 1 ® 1 and S. 
As explained in Section 3, this result allows classifying all stationary states of the 
open dynamics generated by the evolution equations in (7). 
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